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Abstract 

We solve the puzzle of the disagreement between orthogonal polynomials methods and 
mean field calculations for random N x N matrices with a disconnected eigenvalue 
support. We show that the difference does not stem from a Z 2 symmetry breaking, 
but from the discreteness of the number of eigenvalues. This leads to additional terms 
(quasiperiodic in N) which must be added to the naive mean field expressions. Our 
result invalidates the existence of a smooth topological large N expansion and some 
postulated universality properties of correlators. We derive the large N expansion of the 
free energy for the general 2-cut case. From it we rederive by a direct and easy mean- 
field-like method the 2-point correlators and the asymptotic orthogonal polynomials. 
We extend our results to any number of cuts and to non-real potentials. 
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1 Introduction 



Random Matrix Models have been introduced in order to give an approximate statisti- 
cal description of quantum systems involving disorder, chaos, complexity or whatever 
prevents from solving the equations of motion exactly Those models are described by 
a matrix (Hamiltonian, transfer matrix, or scattering matrix) of large size N, which is 
too complicated to be diagonalized exactly, and for which only statistical observations 
of the spectrum are available (see [1, 2] for a review on RMT). 

Most of the quantities of interest and observables are related to the short range 
(in energy scale) behavior of the spectrum (indeed small energies correspond to long 
time evolution, i.e. to equilibrium thermodynamical properties). This is why the short 
range correlation functions are the most studied. 

At first, the simplest models assumed a gaussian weight for the random matrix 
and gave good agreement with observations, provided that the ensemble of matri- 
ces(hermitian, orthogonal, quaternionic.) has the required symmetries (time re- 
versibility,...) [2, 3]. 

It has been observed that the correlation functions of the spectrum possess universal 
properties at sort range, which do not depend on the probability weight, gaussian or 
not. This universality has been proved for a wide range of models by several approaches 
[3, 4, 5, 6, 7], but a very general proof and the exact hypothesis which lead to it are 
still under investigations. 

Moreover, the long range correlation functions appear to share also some universal 
properties, which depend on the probability weight only through a few parameters 
[8, 7]. The most striking example is the 2-point correlation function that we shall 
discuss below. 

In the following we will restrict our attention to the so called Hermitian-One-Matrix- 
Model 5 [2, 9] (hermiticity corresponds to a system with broken time reversibility, for 
instance in the presence of a magnetic field). 

We consider a hermitian matrix M of size N x N with a probability law of the 
form: 

V(M) = e - N ^v(M) 

where V is a polynomial potential bounded from below. 

We wish to study the statistical properties of the eigenvalues (Ai,...,Ajv) of M 
in the large N limit, in particular the density of eigenvalues p(A), and the correlation 
function i?(A,/i), which measures the probability that two of the eigenvalues take the 
values A and /i. 

5 The other ensembles are of course worth considering, but this one is the simplest. 
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Roughly speaking, the eigenvalues tend to occupy a finite interval centered around 
the bottom of the potential well, and in the large N limit, the density of eigenvalues 
p(A) is a continuous function with a compact support. 

The simplest case, where the support is connected, known as the "1-cut case" , has 
been extensively studied [2]. It is then found that the density p(A) is not universal, 
it depends on the details of the potential V, while the connected correlation function 
i? c (A, /i) is universal in the short range regime (|A — fi\ ~ 0(1/N)), but also in the long 
range regime (|A — fi\ ~ 0(1)) once the short range oscillations (of period ~ 0(1/N)) 
have been smoothed out. 

What happens when the potential V possesses several wells, of approximatively the 
same depths? Then the density p(A) has a disconnected support, [a± : bi] U [a 2 , b 2 ] U . . . U 
[a s ,6 s ], each interval [ai,bi] being centered around one well of the potential V. This 
case is known as the multicut (or multiband) case (here s cuts). 

In the multicut case the density is still not universal, whereas the 2-point correla- 
tion function is universal in the short range regime and seems to have some universal 
properties in the long range regime after smoothing: in [10] an explicit form of the 
2-point connected correlation function was given, and is claimed to be universal: in- 
deed, according to the authors of [10] "it depends only on the number of connected 
components of the support and on the position of the endpoints, but not on the po- 
tential". However, more recently several authors [11, 12, 13] have studied the two-cut 
case s = 2: they concentrated on the case of an even potential V (the two cuts are 
thus symmetric [a, b] and [—6, —a]). Using an ansatz for the asymptotic expression 
of orthogonal polynomials in the large N limit, and rederiving the two-point function 
from this ansatz, they observed that the connected correlation function is still uni- 
versal in the short distance regime (which was expected), but more surprisingly, that 
the smoothed connected correlation function in the long range regime depends on the 
parity of iV (N being the size of the matrix). This seems to contradict the former 
result of [10] ! 

In this paper, we will solve this paradox. 

We will show that the semi-classical method of [10] gives the 2-point connected 
correlation function only up to an additional non-universal term, which is already 
present in the free energy, but subdominant at large N in this case. We correct the 
semi-classical argument of [10], and give a simple (and physically appealing) derivation 
of the origin of the additional term, that we compute explicitly. This allows us to 
recover the results of [11, 13] for the symmetric case, and to generalize them to non 
symmetric potentials, without using orthogonal polynomials. 



2 



Using the same semi-classical argument, we are able to derive large N asymptotics 
for the orthogonal polynomials, recovering the results of [11, 13] as well as the general 
s cuts asymptotics which appeared recently in the mathematical literature [14] , and to 
extend these results to the case of complex potentials. 

The effect leading to the new term in the semi-classical calculation is simple enough 
to be explained briefly in this introductory section. 

In [10], the free energy F of the matrix model is derived by a saddle point approx- 
imation. In particular, one has to extremize the action with respect to variations of 
the number rii (i — 1 . . . s) of eigenvalues in each connected part of the support, or in 
other words with respect to the occupation ratio X{ = rii/N: 



dF 

F(x c ) where — — 

ox 



(1.1) 



However, one has here missed the crucial fact that n« = Nxi are not real numbers but 
integers. When Nx c is not an integer, the extremum of F(x) is never reached, and the 
saddle point approximation has to be slightly modified. Roughly speaking the discrete 
sum cannot be approximated by an integral: 

^ e -p(n-iVxc) 2 ^ I dxe -N*g( n f ^ 
n J 

The discrete sum actually depends on how far from an integer Nx c is. For instance, in 
the symmetric case, we have x c — |, and the result depends on the parity of N. 

We will show that (as expected from Eq. (1.2)) in general the result involves elliptic 
theta functions depending on Nx c , thus leading to a quasi-periodic dependence on N. 
This effect is of order N~ 2 for the free energy, but is of order 1 for the computation of 
the orthogonal polynomials and for the correlation functions. It implies in particular 
that there is no regular large N topological expansion (involving only power series 
in N~ 2 ) for the 2-cut matrix model. 

We will find out that the short range correlation function is universal, while the 
long range smoothed correlation depends on N quasi-periodically. 

The paper is divided as follows: in section 2 we introduce the method and notations 
for the 2-cut model, and we compute the free energy. In section 3 we derive the 2-point 
correlation function, and we recover the expression of [11, 13, 14] in the symmetric 
case. In section 4, we give an asymptotic expression for the orthogonal polynomials, 
which we use to rederive the universal short range properties of the spectrum, as well 
as the smoothed long range 2-point correlation function. 

The generalizations to a complex potential or to an arbitrary number of cuts are 
presented in Appendix B and C. Appendix A is a summary of some relationships 
between elliptical functions in case the reader is not familiar with them. 
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2 The free energy 
2.1 Basics 

We start from the standard Hermitian matrix model defined by the partition function 

Z[V-N] = J d N [M] e - NtlV M (2.1) 

where N is the dimension of the matrix M, V is an analytic - in general polynomial - 
and for the moment real - function, and djv[M] is the standard U (N) invariant measure 
over Hermitian matrices 

N 

d N [M] = \[dM ri J] 2dRe(M ij )dlm(M ij ) (2.2) 

i=l l<i<j<N 

Integrating out the " angular part" of M, Z can be rewritten as an integral over the N 
eigenvalues \ 1: . . . , \ N oi M [2] 

Z[V-N\ = C N Z[V;N] (2.3) 

/N _ . N 

Y[d\ k e~ N ' EkV(Xk) JJ(A fe -A;) 2 = Y[d\ k e- N2s ^-^ (2.4) 

k=l k<l J k=l 

with the measure factor 



C N = Vol 
and with the action S(Xk): 



XJ(N) 



v(i) N x e N 



J_ fr (^y^_ , 25) 



1 N 1 

5(Ai,...,A^) = -X;V(A*)-2^ £ ln (^-^) 2 (2-6) 

fc=l l<k&<N 

In the simplest "one cut" case, corresponding in particular to a concave potential, it is 
known [9] that the free energy F defined as 6 

Z[V;N] = ( 2 -Pi 2 e"^ (2.7) 



has a topological large N expansion 

F = N 2 F + F l + N~ 2 F 2 + . . . (2.8) 
6 Z is normalized here so that F is zero for the Gaussian model V = \M 2 
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obtained for instance by re-organizing the perturbative expansion according to the 
topology of the Feynman diagrams. The large N limit (planar limit) can be described 
by a "master field" configuration where the eigenvalues are described by a continuous 
density p(A) with a connected compact support C = [a, b], with the constraints 



dAp(A) = 1 and p(A) > if A G C (2.9) 

and the action 2.6 becomes 



c 



S[p] — f d\V(\)p(\) - I d\dp p{\)p{p)\n\X- p\ (2.10) 
Jc Jcxc 

The leading term of the free energy F can be obtained by the saddle point method: 
the effective action S[p\ is extremized for a continuous distribution p c and we have 
simply 

F = S[ Pc ] (2.11) 

(up to an additive - potential independent - constant). To compute p c we include the 
constraint 2.9 in the effective action by a Lagrange multiplier T 

S[p] = S[p] + T(l-^P) (2.12) 

The saddle point equation for p reads: 

——- = V(X) - 2 / dpp(p) ln(|A - ^1) - T = VAGC (2.13) 
<9p(A) J c 

which simply means that the real part of the effective potential 

Kfr(A) = V(X) - 2 J dpp(p) In (A - p) (2.14) 

is constant on the e.v. support C, and equal to T. The derivative of 2.13 w.r.t. A gives 
the well-known equation 

Re(cu (A)) = Idpp(p)—!— = V'(X)/2 VAgC (2.15) 
J c A — p 

where u is the large N resolvent 



A-M 



= /d^fW (2.16) 

JC A — p 



Finally let us recall that in the one-cut case C = [a, b], if the potential V is a polynomial 
of degree P, uj is of the form 

= _ Mwfm with = _ _ 



where M(A) is a polynomial with degree P — 2. a, 6 and M are entirely determined by 
the constraint that 



uo Q (\) = \- 1 + 0{\- 2 ) for A^oo (2.18) 
The e.v. density is given by the discontinuity of cu 

p(a) = I mx + i 0+) - lo(x - io +)] = m^m (2 . 19) 

2.2 The 2-cut case 
2.2.1 Mean field: 

If the potential V is real but has more than one minimum, the large N limit may be 
described by an e.v. distribution on several disconnected intervals. For simplicity we 
shall first consider the case where there are two intervals 

C = Ci U C 2 , C 1 = [a,b] , C 2 = [c,d\ , a < b < c < d (2.20) 

In this case, as we shall see, there is no topological large iV expansion, even for 
the free energy F. As shown in [15, 16], to describe the large iV limit, we have to 
consider as an additional variable the "average" proportion of eigenvalues Xi = ni/N 
and x 2 = n 2 /N in each interval C\ and C 2 , and introduce the associated Lagrange 
multipliers Ti and T 2 for the constraints 



x a = / p(A)dA , a = 1,2 (2.21) 
Jc a 

The effective action 2.12 now reads, with 

x = x x (2.22) 

S[p; x] = S[p] + J2 T «( X «- I P( X ) dA ) . x 1 + x 2 = l (2.23) 

a=l 

with S[p] given by 2.10 as before. The saddle point equation w.r.t. p(A) gives as before 
the equation 2.13, which implies that the effective potential defined by 2.14 is constant 
on each interval 

V eS (X) = T a when A e C a (2.24) 

but the corresponding e.v. density p c (X) and the effective action S c still depend explic- 
itly of the e.v. proportion x, since we have 

S c [x] = \ Pc (X)V(X)dX + ^r a x a \ (2.25) 
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The saddle point equation w.r.t. x implies the equality of the effective potentials for 
each interval 

= r\ - r 2 = o (2.26) 

This fixes the value of x, and it is known that with this last equation the e.v. density 
p c is uniquely determined in the 2-cut case. The large N free energy is then given 
simply by 

F = S[p c ;x c ] = S c [x c ] (2.27) 

For an explicit polynomial potential V of degree P, and for fixed x, the 2-cut mean 
field solution for the resolvent is 

oj (X,x) = - M(A) V^^ W ith a(X) = (A - a) (A - 6) (A - c)(A - d) 

(2.28) 

and M(A) a polynomial with degree P — 3. The e.v. density p(X,x) is still given by 
the discontinuity of ujq. The coefficients of M and the 4 end-points a, b, c, d are entirely 
determined by the constraint that ouo(X) ~ A -1 when A — > oo and by the fact that x 
must be given by 

x = J\(\,x)d\ = ^ J\m{\x)\^\^\<\\ (2.29) 
Finally the equation 2.26 which fixes x = x c reads 
= V cS (b)-V cS (c) = f dX (2uj (X,x)-V'(X)) = - [ dXM(X,x)y/\a(X)\ (2.30) 

Jb Jb 

2.2.2 Discreteness of number of e.v.'s: 

This is sufficient if one is interested in the leading term in the large iV limit (planar 
approximation). However, in order to understand the structure of the subdominant 
terms of the large iV expansion, it turns out that we cannot neglect the fact that the 
number of e.v. n a = Nx a in each interval C a must be an integer. 

Let us consider the simple case where the potential V has two separate minima z\ 
and z 2 . What has to be done is first to fix the number of eigenvalues n\ — n (resp. 
n 2 = N — n) in the vicinity of z\ (resp. z 2 ) in the partition function 2.5 by writing 

Z[V;N] = j: "l n) , Z[V-n,N-n\ (2.31) 

where 

/E _ r+oo 

00 i<n J E j>n k<l 
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with E a "frontier" b < E < c between the two semi-classical cuts [a, b] and [c, d]. We 
now claim that each term of this discrete sum has a well defined large N topological 
expansion. Indeed, we can rewrite 2.32 as a matrix integral over two separate matrices: 
a rii x n \ matrix Mi with the n± — n e.v. < E and a n 2 x n 2 matrix M 2 with the 
n 2 = N — n e.v. > E 1 , as 

= — \ /d ni [Mx] /■ dn jM 2 ]e-^^^ Ml »-^^^ M2 )) +2 ^( In ( M ^ Id+Id ® M2 )) 

^n^N-n J J 

(2.33) 

This last matrix integral has a topological large iV expansion of the form 2.8 in the 
't Hooft limit N — > oo, x = n/N fixed, obtained by doing a classical perturbative 
expansion around the smallest minimum z\ of V for Mi and around the largest mini- 
mum z 2 of V for M 2 , and by re-organizing the perturbative expansion according to the 
topology of the Feynman diagrams. Taking into account carefully the measure factors 
C n and Cjv-n, an d using their large N asymptotics 

C N = ^ f^j 2 e % N2 (2n)- N N& est (1 + C(^ 1 )) when iV ^ oo (2.34) 
(easily derived from Stirling formula), we obtain that 

N 

N -±J2e~ F ^ N ^ (2.35) 

n=0 

where each F[V; N,x] has a regular large iV asymptotic expansion of the form 

oo 

F[V;N,x] = N 2 - 2h F h [V,x] where x = n/N (2.36) 

h=0 

with each F h [V, x] a regular function of x — n/N. In particular, the leading large N 
term is given (up to an additive - V and x independent - constant) by the classical 
effective action 2.27 

F = S[p c ;x c ] = S c [x c ] (2.37) 

Finally, let us stress that although this decomposition depends on the arbitrary param- 
eter E, since E is in the interval ]b,c[ where the density of eigenvalues is exponentially 
small with N, the integral 2.32 depends on E only through exponentially small terms 
of order e ~ cst Ar , which are "non-perturbative" in the topological expansion Eq. (2.36). 

2.2.3 Beyond mean-field: 

We can now easily calculate the sub leading terms of order 0(N~ 2 ) for the full partition 
function. In the large N limit we can approximate the sum 2.36 by 

N 

Z[V;N] oc Y^^ n2Fo[V ' M - Fi[V '' x]+ - (2.38) 

n=0 
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If x c denotes the saddle point of F [x] given by 2.26, the sum is dominated by the n's 
such that 

\n-Nx c \ = 0(1) (2.39) 
Thus we can still use a quadratic approximation for F [x] 

Z[V;N] OC e -{ N2F o[ V M+Fi[VM + -) ^ e -{n-Nx c fF{][V;x c ]/2 ^AO) 

n 

where Fq = d 2 F /dx 2 and where the • • • represent terms of order 0(N~ 2 ). The last 
sum over n gives simply an elliptic Jacobi theta function 9 3 

j2 e -(n-N Xc )*F>>iv,x c ]/2 = ^nF^[V; x c ]y 1/2 9 3 (Nx c \r) (2.41) 

n 

with modular parameter r given by 

r = ^ j (2-42) 
and where the theta function is defined as 

3 (z\t) = 9 3 {z) = J2q n2 e 2iwnz with q = e mT (2.43) 

It obeys the periodicity relations 

6 3 (z+l) = e 3 (z) , e 3 (z + r) = e-^ 2z+ ^e 3 (z) (2.44) 

(For details on elliptic functions see e.g. Refs [24, 25, 23]). Eventually we have for the 
free energy 

F[V;N] = N 2 F [V,x c } 

- In (6 3 (Nx c )) + F 1 (V; x c ) + \ In (2nF^'[V; x c }) (2.45) 
+ 0{N~ 2 ) 



where Fi is the torus contribution in the topological expansion of 2.33. The next terms 
of this expansion can be calculated along the same line. 

Let us stress that this is not a topological expansion, since the second term In (#3(iVx c )), 
seemingly 0(1) and contributing at the torus order, is not regular in N. Indeed, it is 
periodic in x c with period 1/N. When computing some observables or quantities of the 
matrix model, one must take derivatives of F w.r.t. some parameters of the potential 
V . Since the saddle point x c depends implicitly on V, every derivative will give a factor 
N, and this term may become of the same order than the first term N 2 F [x c ] given by 
the planar limit. Note that the last two terms depend on x and not on Nx, and they 
will remain subdominant once we take derivatives of F. 
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u(c)= T/2 



u(b)=l/2+ t/2 



a b 



c d 



u(d)=0 




u(a)=l/2 



Figure 1: the upper half-plane is mapped onto a rectangle (1/2, t/2) 



2.3 The modular parameter 

Finally, we can express simply the modular parameter r defined by Eq. (2.42) in term 
of the end-points a, b, c, d of the support of e.v. For this purpose, we introduce the 
function a 

(7(A) = (A - a){\ - b){\ - c)(A - d) (2.46) 



and the function u 

where K is 

pc dz 



dz 



(2.47) 



K = 



K[m] with m = ( , d ~^ C ~ b l (2-48) 



/ 6 V^(^ y/(c-a){d-b) (d-b)(c-a) 
K[m] is the complete elliptic integral of the first kind. Similarly we define 
dz 2 



K' = 



a V\ a ( z )\ y/(c-a)(d-b) 



K[m'} with m' = 1 - m (2.49) 



We shall show that the modular parameter r of Eq. (2.42) coincides with the stan- 
dard modular parameter of the torus associated to the mapping u, i.e. of the elliptic 
curve y 2 = <j(z). Indeed, r is simply given by 



.IT K[l-m] 
T ~ l ~K ~ 1 K[m] 



(2.50) 



So we have 
u(d) = 



u(b) 



1 + r 



u(c) 



- , u(oo) = Moo 

(2.51) 

and u maps the upper half A-plane onto the half-periods rectangle (1/2, t/2) and the 
double-sheeted complex A-plane onto the period rectangle (l,r). 
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To show Eq. (2.50), we use the fact that in the two-cut case, if we fix x (the e.v. 
ratio in the first cut) the semiclassical e.v. density (extrema of the effective action 
S ) is now a function p(A, x) of A and x, and the end-points a, b, c, d depend on x. 
Therefore the large N resolvent uj is of the form 



aiM .m_M (2 . 52) 

with M(A, x) a polynomial with degree P — 3 in A (P being the degree of V), entirely 
fixed by the constraints 2.18 and 2.29. Therefore the partial derivative of c^o(A, x) w.r.t. 
x is necessarily of the form 

du (\x) = C 
dx y^(A) 

with C = C(X,x) a priori a polynomial in A. Since 2.18 still holds independently of x 
we must have 

duJ i X ' x) = 0(X- 2 ) for A - oo (2.54) 
ox 

which implies that C(X,x) is of degree in A, i.e. is a constant (depending only on x) 

C = C(x) (2.55) 
This constant can be easily determined by using that 

x — p(X,x)dX = / — uj (X,x) (2.56) 
J a Jc> 2ivr 

with C a clockwise contour encircling the interval [a, b] . Therefore we have 

£ = 1 = / = => C = -JL ( 2 .57) 

with K' the half-period defined in Eq. (2.49). Now we use Eq. (2.26), Eq. (2.27) and 
the definition of the effective potential V e g of Eq. (2.14) to write the derivative of the 
free energy w.r.t. a; as 



dF 



V cS (b) - V cS (c) = f dX (2uj (X) - V'(X)) (2.58) 
Jb 



dx 

Now we take the derivative w.r.t. x of this equation and obtain 

= = 2 TdA^l = 2 TdA^ = -2CK = 

dx 2 J b dx J b \fo\X) K> 

Using Eq. (2.42) we thus obtain the result 2.50. 



(2.59) 
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3 2-points correlation function 



3.1 The basic formula 

As a first application we compute the large N smoothed connected two-point correlation 
function (first obtained by [11, 13, 14]), defined as 



cu c (A,/i) = (Tr 



1 



X-M 



Tr 



p-M 



- (Tr 



X-M 



)(Tr 



jx-M 



Adding source terms to the potential of the form 

V ex = V(z) - e xj ^ , V ex ,^(z) = V{z) - e x 1 



we have 



u c (X,n) 



X-z 



li — z 



d d 



iV 2 dexde^ 
and for the resolvent (one-point function) 



eA=en=0 



u(X) 



1 

N 



Tr 



1 



X-M 



_L_ _8_ 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



e x =0 



If the e\s are small and the A's not too close to the cuts, the mean-field solution is still 
a two-cut e.v. distribution, with x c = x c (e) an explicit function of the A's. So from 
Eq. (2.45) for the free energy we have for the two-point function 



u c (X,fi) = 



d d 
de x <9e M 



Fo[Ve^] + ^^[ln(e 3 (Nx c ))}" 



+ 0(N- 1 ) (3.5) 



The first term in the r.h.s. of 3.5 is the mean-field contribution already calculated in [17, 
10], the second term involving a second derivative of an elliptic function, characterizes 
the multi-cut solution. 



3.2 The mean-field contribution 

For completeness let us first rederive the mean field contribution of [17, 10]. Taking 
the derivative with respect to e A we obtain the mean-field resolvent for the potential 

d 



de x 



F [V e 



^o(A;KJ 



e A =0 



which must be of the form 



^o(^;KJ = 



V'(z) 



(z - fif 



M(z)yft(zj ' 

(Z - »y 



(3.6) 



(3.7) 
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with M(z) a polynomial of degree P — 1 (here both the coefficients of M and of a 
depend on /j, and e M ). In addition to the P — 2 constraints (coming from 2.18 and 2.30) 
u)q(z; V 6fi ) must be regular at z — /i. This determines entirely M. Taking the derivative 
w.r.t. e M and using the symmetry A <-> /i we get 



wg(A,ji) = 



de x de^ 



F n = -- 



2 (A - /j,) 2 



1 + 



Q(A,/x) 



y/a(X)y/a(ti) 



(3.8) 



with Q(X,/jl) a symmetric polynomial in A and in /i. The constraints on are: (i) 
= 0(\~ 2 ) as A — > oo which implies that Q is of degree at most 2; (ii) tug is regular 
at A = /i which implies that Q(X, /i) = — er((A + /i)/2) + 0{{\ — /i) 2 ); (hi) finally the 
equality of the effective potential on the two cuts implies that 



(Vq(\, /i) d/i = 



(3.9) 



Conditions (i) and (ii) fix uniquely Q 



Q(A,/i) = - 



(A — a)(/i — b)(/i — c)(A — c?) 
+(// - a) (A - 6) (A - c)(/i - d) 



+ 5(A-a«) ; 



(3.10) 



up to a constant 5 fixed by condition (iii), which is found to be 

(c — b)(d — a) 



S = -2 {c - a){d - h) K[ri 



with 



m 



a){d — b) 



(3.11) 



and where K[m] and E[m] are the standard elliptic integrals of the first and second 
kind. 



3.3 The non-regular contribution 

In order to compute the second contribution to 3.5, we simply need |ff . Since x c is 
fixed by the constraint ^ = we can write 



8x r 



d 2 F n 



de\ dxde\ 
Using the results of subsection 2.3 we have 
2ttK 



' d 2 F 
dx 2 



d 2 F 



dx 2 K' 
So we have eventually 



d 2 F d 7T 1 

dxde x dx K' y/a(X) 



dXr 



de x 2K v fc(Xj 



(3.12) 



(3.13) 



(3.14) 
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and the second non-regular term is 
0j '' ° Xn [\n(0 3 (Nx c \r))r 



2K^{X) 2K^f^) 



[\n(e 3 (Nx c \r))}" (3.15) 



with K defined by Eq. (2.48). Using standard relations on elliptic functions, this can 
be rewritten as 



(c — a){d — b) 



4 V ^(A) V ^) 



E[m\ 
' K[m] 



dn 2 (iVx c + i) 



(3.16) 



with 



dn(w) = dn(2K[m]u\m) (3.17) 

where dn(w|m) is the Jacobi elliptic function dn. Its periods are 2K[m] and 4ilf[ra'], 
and dn 2 (z) has periods 1 and r. 

3.4 The final result 

Combining 3.8, 3.10, 3.11 and 3.16 we obtain the final result for the 2-point correlation 
function 



w c (A,/i) = 



1 



4(A-/i) 2 
(c — a)(d — b) 



(A — a) (A — b)(fi — c)(/i — d) 
(fj.-a){fj.-b)(\-c){\-d) 



+ (A <-> n) 



sn 2 (iVx c + ±) 



(3.18) 



4 V ^(A) V ^) 

We have used the relation dn 2 (w) = 1 — msn 2 (w), where similarly to 3.17 we note 

sn(u) = sn(2K[m\u\m) (3.19) 

Surprisingly, the ratio E[m}/ K[m\ characteristic of the mean-field solution of [17, 10] 
has disappeared. 

The smoothed 2-point connected density correlator p c (A,/x), defined as 

p c (A, ijl) = ( Tr [5(X - M)\ Tr [5(jm - M)]> - ( Tr [5(X — M)])( Tr [5(jjl - M)]> (3.20) 

can be obtained easily from the discontinuity of u c (X,fi). One obtains in the large N 
limit, if A and \i are on the support of e.v. 



p c (A,/x) 



4vr 2 



(A — a) (A — b) (/i — c)(fx — d) 



(A - A*) 2 



(^-a)(^-6)(A-c)(A-d) 



c — a)(d — b) 9 ., r . 

sn 2 (iVa; c + ±) 



+ A <-> /x 
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(3.21) 



s x = l if Ae[a,b] , -1 if Ae[c,d] , (3.22) 
and zero otherwise. 

The new non- regular term sn 2 (iVx c + ±) is an even periodic function of Nx c with 
period 1 which varies between and 1. Therefore, as iV varies, depending on the 
rationality or the irrationality of x c , the two-point function will be varying with N in 
a periodic or quasiperiodic way. 

3.5 The symmetric case 

It is now very easy to recover the results of [11, 13] for a symmetric potential. Indeed, 
if the potential V is symmetric, the two cuts are also symmetric 

a = -d , b = -c (3.23) 

and we have automatically 

x c = \ (3.24) 

so that 

2/ at x f sn 2 (i) = 1 if N is even 

Sn2 (^c+|) = { sn 2$ = 0if7Visodd ^ 



Eq. (3.18) and Eq. (3.22) become 



2 



(-1)" a± 
2 v^Alv^) 



/(A, /*) = A n i ( (a2 "y ( 'l" V) -(-If ab] (3.27) 



2vr 2 ^^A^^^/jf V (A - /i) 2 
with (7(A) = (A 2 -a 2 )(A 2 -6 2 ). 

3.6 The two-point function as an elliptic function 

It is interesting to consider the two-point correlator in terms of the elliptic coordinates 
defined by Eq. (2.47) 

u = u(A) , v = u(fi) (3.28) 

Let us thus consider 

cu c (u,v) = = 2K v ^(A)2^ v ^)^ c (A,/i) (3.29) 

It is easy to see (from the properties of uj c ) that Cu c (u,v) satisfies: 

1. Q c (u,v) is a doubly periodic function of u (and of v) with periods 1 and r; 
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2. uj c {u,v) is regular at u and v = u(a), u(b), u(c), u(d) and Uoo] 

3. uj c {u,v) is regular when u — v, but has a double pole at u — —v (corresponding 
to the double pole of ou c (\, //) when A = \i but with A in the first sheet and \i in 
the second sheet), with residue 1. 

This implies that u c (u, v) is a Weierstrass elliptic function 

uj c (u, v) = p(u + v\t) + constant (3.30) 

where the constant depends on Nx c (p has periods 1 and r). Using classical identities 
between the Weirstrass p function and the Jacobi elliptic functions, it can be easily 
calculated. We find the remarkably simple result 

co c (u,v) = p{u + v\r) - p(Nx c + ^\t) (3.31) 

or equivalently 

lu c (u,v) = - pn(0i(u + u|r))]" + [In (9 3 (Nx c \r))}" (3.32) 

4 The orthogonal polynomials 

Let us briefly recall some basic facts about the well-known method of orthogonal poly- 
nomials [18], which is a powerful tool for studying the spectral properties of random 
matrices [2]. Asymptotic expressions for the orthogonal polynomials have been ob- 
tained recently [14] in the mathematical literature, by solving a Rieman-Hilbert prob- 
lem. Here we will derive them from the free energy directly. 
Consider the partition function (2.4): 

Z = I dAL.-dAtfe-^W n^-A,) 2 (4.1) 

J i<j 

The last term is a Vandermonde determinant [2]: 

T[(\ t - A,) = det ((A,)^ 1 ) = det {V^K)) (4.2) 

Kj 

where the last equality is obtained by linearly mixing columns of the determinant, and 
holds for arbitrary monic polynomials V n (X) with leading coefficient V n (X) = \ n + . . .. 

The method of orthogonal polynomials consists in choosing a family of polynomials 
suitable for the computation of (4.1), namely, the family of polynomials orthogonal 
with respect to the weight exp — NV(X): 

J dA V n {\) V m (X) e' NV ^ = h n 5 nm (4.3) 
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With this particular choice of polynomials, the integral (4.1) is merely: 

JV-l 

Z = N\\[h n (4.4) 

n=0 

and the joint probability density of all the eigenvalues takes the form of a Slater de- 
terminant: 

R N (X U ...,\ N ) = -L jet N [Vn-i(A,)] (4.5) 

\l<i<JV / 

where the wave functions ip n (\) = -j=V n (X)e~% v W are orthonormal. 
4.1 The Kernel K(\,/j) 

The square of a determinant can be rewritten as the determinant of a product: 



det (ip n -i(\)) = det 

n,i I l<i,j<N 



2 



'JV-l 



5^n(Ai)^„(Aj) 



n=0 



we are thus led to introduce the kernel K(X,/i) [19]: 



JV-l 



n=0 



In terms of which the joint density of eigenvalues is now a determinant 

N N 

R N (X U . . . , X N ) = — det [K(X h Xj)] (4.7) 

The orthonormality properties of the polynomials imply the projection relations 

J dA K(X, X) = 1 and J dA A)X(A, i/) = u) (4.8) 

which make any partial integration of (4.7) easy to perform (theorem of Dyson [20]). 
In particular, the integration over N — 1 eigenvalues gives the density of eigenvalues 

p(Ai) = J dA 2 . . . dA ArJ Rj V (Ai, ...,Xn) = K(X h Ai) 
and the integration over N — 2 eigenvalues gives the correlation function: 



# 2 (Ai,A 2 ) = /dA 3 ...dA A r J R iV (Ai,...,A 



JV 



JV 



(if(Ai, Ax)K(A 2 , A 2 ) - tf(Ai, A 2 )K(A 2 , Ax)) 
In short: 

p(A) = X(A, A) , p(A, p) = (K(X, X)K(fi, y) - X(A, /i) 2 ) (4.9) 
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In addition, the Darboux-Christoffel theorem [24, 21], asserts that 

Nh,N—\ A — /i 

which means that we need to evaluate V n only for n = N and n — N — 1. 

Thus, we shall now aim at finding asymptotic expressions for the orthogonal poly- 
nomials V n (\), and the kernel K(\,/j,) in the large N limit, and n close to N. This 
has been done in the 1-cut case [7] and in the symmetric 2-cut case [11, 14, 12]. Here 
we will generalize it to the non-symmetric case, with the method used in [21]. 

4.2 WKB approximation for the orthogonal polynomials 

The orthogonal polynomials have the following integral representation (see Appendix 
1 of [21] or [18]): 

v(x) JdM nxn det(A-M)e-^W 
' n[) ~ JdM nxn e- N *nM) 

where the integral is restricted to hermitian matrices of size n x n. 

Thus the orthogonal polynomial is given by the ratio of two matrix integrals of the 
same type as the partition function 2.1: 

Z[V + 5V 1 + 5V 2 ;n] _ e -^W 1+ ^ ; n] 
' n[) ~ Z\V + bV x -n\ e -F[v + s Vl ;n] ^- LZ ) 

where 

N — n 1 
5V 1 (z) = V(z) and 5V 2 = In (z - A) (4.13) 

71/ 71/ 

We have seen in the previous section (eq.2.45) that 

F[V;n] = n 2 F [V;x c ] - \n6 3 (nx c [V}) + ... (4.14) 

We will use the fact that under a variation SV of the potential, the variation of F is 
[21, 8]: 

5F = — I u(z)SV{z)dz (4.15) 
2m J 

where the anti-clockwise contour encloses the support of the density of eigenvalues, 
and u(z) is the resolvent (eq 3.4 and 2.28): 

U(Z)= 1 -(V'(Z)-M(Z) V ^(Z) 

It is convenient to introduce two sources t\ and t 2 for the variations SVi and 5V 2 of 
the potential, and consider a generalized potential V(z): 

V(z) = V(z) + hSV^z) + t 2 5V 2 (z) = V(z) + hV(z) + t 2 In |A - z\ 
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Since t\ = and t 2 = — \ are both small of order 0(N 1 ), we will expand F in 
Taylor's series: 

t 2 t 2 
F[V(z); n] = F[V(z); n] + t&F + t 2 d 2 F + -±dnF + t 1 t 2 d 12 F + -^d 22 F + ... 

all the derivatives being taken at the point t± — t 2 — 0. 
This will give 

V n {\) ~ e^e^-^^e-i^ + ( N & x ) x 

v 7 6> 3 (nx + (TV — n)9ix) 

Now, let us compute the derivatives of F and x = x c with respect to t\ and t 2 - 
The method proceeds similarly to section 3.1. 

4.2.1 Derivatives of F with respect to t\ and t 2 

using (4.15) with (4.13): 

dF Q 



/u(z) In (z — \)dz 



dt 2 2i% 

After integration by part, the pole in [z — A) picks a residue, and the result is a primitive 
of u(X): 

dF 



u(z)dz (4.17) 



The lower bound of integration Aq is to be chosen such that e nd2F ° ~a-+oo A™, i.e. 



f°° 1 

lnA = / (u(z)--)dz (4.18) 

In order to compute the second derivatives 8\ 2 Fq and d 22 Fo, we will need to differ- 
entiate cu(z) with respect to ti and t 2 . 

4.2.2 Derivatives of iv(z) with respect to t\ and t 2 

The resolvent a; (2) computed for the potential V(z) takes the form: 

u(z) = ±(y(z)-M(z)yfiW) (4.19) 

where M(z) is analytic. Notice that when V'(z) has a pole in z = A, M(z) may have 
a pole too. cu(z) obeys a linear equation: 

u{z + iO) + u{z - iO) = V(z) for z e [a, b] U [c, d] (4.20) 

Thus its derivatives obey linear equations as well: 

9icj(^ + iO) + diw(z - iO) = 5V[{z) = V'(z) (4.21) 
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d 2 uo(z + iO) + 9 2 cj(^ - iO) = £V2(z) = — *— (4.22) 

2 — A 

• (9cji: The solution of (4.21) is: 

<M*) = - 4= (4-23) 

where /(z) is analytic in z. The boundary conditions 2.18 imply that f(z) ~ z when 
2; — > 00 and / has no pole, thus f(z) is a polynomial of degree 1: 

<4 ' 24) 

is determined as a function of a, 6, c, d by the derivative of 2.30 with respect to ti. 

dz l^^ = (4.25) 

It can be checked that in term of elliptic theta functions we have (see Appendix A, or 
[23]): 

z-z _ d ln ^Kz)+M 

VM-Z) ^ ^1(^(2) - Moo) 

and thus: 

^M*) / \ d 9i(u(z) +M00) n _v 

= ufz) - — In - ; ; { f (4.27) 

dh y ' dz 0i(u(z) - Moo) 

• cta^: Note that the t 2 source-term is the primitive of the eA source-term of 3.2, and 
that 

d du(z) du(z) 1 <9 2 F 

Tz~dto~ = ~d^ = -n*d^dT x = ^ A) (428) 
so, du/dt 2 has already been computed in 3.8. The second derivative d 22 F corresponds 
to z — A. 

d 22 F = In v^A) + 2 In (0i(u(A) - ««,)) (4.29) 

4.2.3 Derivatives of x 

Recall that 

using 4.27 we get: 



f p(A)dA = -L T M(*) xAMdz (4.30) 



— = x + — / -i===dz = x + 2 UoD 4.31 

Similarly, from 4.28 and 3.18 or 3.31 (or less tediously, taking the primitive of 3.14), 
we get: 

(FT 

^- = -u(\)+ Uoo (4.32) 
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4.3 Final result 
4.3.1 Case A [a, b] U [c, d] 

Eventually, inserting 4.17,4.27,4.29,4.31, 4.32 into 4.16 we get: 



A$[a,6]U[c,<i] 



(4.33) 



where 



p n («) = C n 



9 3 (Nx + 2(7V - n)^ +u- m 00 )6'i(2m 00 ) f O^u + u 



9 3 (Nx + 2(N - n^O^u - u c 



QAu - u c 



n-N 



(4.34) 



C n is a normalization such that V n ~ A™ for A — > oo. 

1 W) 



2fC4 



n-JV+1 



with 



2tf0 1 (2u oo ) 

3 (O) 



and 



^ |tZ- a - c + 61. . 

4 ^ 3 (2moo) 



(4.35) 



Note that 4.33 is unchanged under u — > u + 1 and u —> u + t. Indeed, a shift 
■u — > u + r amounts to a nontrivial circle around the cut [c, d] . Thus J* uj is shifted by 
— 2«7r f^p = — 2«7r(l — x), and e N receives a phase e 2Mr7V:c . In the same time, the 9 
functions receive phase factors: 9(v + r) = 0(i;)e -2 " r ( r ' +T / 2 ). One can easily check that 
the total phase shift is 0. 

4.3.2 Case A e [a, b] U [c, d] 

Expression 4.34 has been derived by a saddle point approximation of 4.12 when A does 
not belong to [a, b] U [c, d}. When A lies on the cut [a, b] U [c, d], 4.12 actually has two 
saddle points, contributing to the same order. They correspond to the two determi- 
nations of the square root ±a/o"(A). The asymptotic expression for the orthogonal 
polynomial is then given by a sum of two terms: 

VJX) = CVti\pJu)e- iN ^) + i Pn (-u)e iN ^] e^ y(A) (4.36) 

\e[a,b]L>[c,d\ J 



where ((X) = p(z)dz and 



C 



d N e 2 



f V (d)-Nf™(u(z)-l)dz 



(4.37) 



To summarize: When A [a, b] U [c,d], the wave function ijj n (X) = V n (\)e~ NV ^ 2 
decays exponentially, and within the support [a, b] U [c, d], it oscillates at a frequency 
of order N. 
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Figure 2: Typical behavior of the wave function 



4.3.3 Check of the orthogonality 

For completeness, let us check that the functions (4.33) are indeed orthogonal (at 
leading order in N' 1 ). Let us compute the integral: 



/oo 
dXV n (X)V m (^ NV{X) 
-oo 



The contributions of the integral along ] — oo, a] U [b, c] U [d, oo[ are exponentially small 
and do not contribute at leading order. 

Along [a, b] U [c, d], we use expression 4.36, and get a sum of four terms: 



C 2 / dA u'(A) 



iPn{u)p m {-u) 



+ip n (-u)p m (u) 



(4.38) 



+Pn(«)p m («)e- 2lA ^ A ) -p n (- M )p m (- M )e 2 ^ A ) 

Since the two last terms have fast oscillations of frequency N, they are suppressed as 
0(1/N). 

The leading contribution is thus given by the two first terms of 4.38, which can be 
rewritten as integrals in the u plane along the contour depicted on fig. 3. a: 



/ 



-NV 



= IC. 



du 



6 3 (x n + U - U^O^Xm - U - Moo) f &l(u + Wqo) 



9?,{Xn)0z{x m )9i{u - UozjOiiu + Moo) \0l(u - U c 

+ (u — > — u) 



(4.39) 



where x„ and c„ m are short notations for: 



x n = Nx + 2(N - n)Moo and c nrn = C 2 C n C m 6' 2 (2Moo) 



(4.40) 



If n > m we may deform the contour to a circle around the point —u^ (fig.3.b), 
and the integral vanishes since there is no pole, while if m > n we deform the contour 
to a circle around +W00 (fig-3.c). Therefore, the integral vanishes for n 7^ m. 
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Figure 3: Deformation of the contour integral 



When n = m, the integral picks a residue: 

J d\V n V m e~ NV ^ = hj nm (4.41) 
with (C, c nm , A, defined in 4.35, 4.40, 4.37): 

h - C — 9 ^ X n+l) _ _ 4 C 2 Aafon+l) A 2(n~N+l/2) U 42 ) 

nn - Cm 6 1 (2a oo )ff 1 (0) 9 3 (x n ) - 47r ° 9 3 (x n ) A ^ 
4.3.4 Recurrence equation 

It is well known that the orthogonal polynomials satisfy a recurrence equation of the 
form [2, 18]: 

XV n (X) = V n+1 {\) + p n V n (X) + a„P„_i(A) (4.43) 
Here, we find that (divide 4.43 by V n , and match the poles on both sides): 

_ h n _ 2 6> 3 (:r n+ i)6>3(:r ra _i) 

which can be rewritten more compactly as 

Oi n = 4 - a) - (c - 6)) 2 + 4(d - a)(c - b) cn 2 {x n + 1/2) ) (4.45) 
lb 
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And by taking u = in 4.43 we get f3 n : 
P n -d = A 



6z{Xn+Z/2)Qz{Xn) 9 3 (x n+1 )6 3 (x n _ 1/2 ) 



_6 3 (x n+ i)6 3 (x n+1/2 ) 3 (x n )6 3 (x n+1/2 ) 
which can be rewritten more compactly as: 

a + d+(c-b) f _ n d-b 

cn 2 (a; n -Woo+l/2) 



(4.46) 



A. = ^ {c ~ b) ^b^ ^ (447) 



The sequences a n and f3 n are thus quasi-periodic in n. It is interesting to recall that 
the behavior of these coefficients has been extensively studied (mainly by numerical 
methods) by several authors in the early 90's [22]. The general conclusion was that 
in the multi-cut case the general behavior of the recursion coefficients was "chaotic" 
in n (and regular or quasi-periodic only in some special cases). It is clear from our 
expressions that in the two-cut case the behavior is always periodic or quasi-periodic 
and never chaotic (in the mathematical sense). This is in fact true even if the number 
of cuts is larger than 2 (see appendix C). 

In the symmetric case, x — \ and = |, we have x n = n/2 modi, so that we 
recover f3 n = and a n = \{a — (— l) n b) 2 . 

In the general case, a n and j3 n vary along a periodic curve, between two extrema, 
given by: 

(d-a-(c-b)) 2 ^ ^ ( d -a + (c-b)) 2 

< a n < 

16 ~ ~ 16 

d + a c — b d + a c — b 

— Pn — ^ I" 



2 2 ~ ' 2 2 

Similarly to the one-cut case, one may relate a n to square width of the distribution 

of eigenvalues, and f3 n to the center of the distribution. 
4.4 The kernel K(X,/i) 

We can now evaluate the kernel K(X, /i) according to (4.10). Let us note u = u(X) and 
v = u(n) and we assume A, /i e [a, b] U [c, d\. 

K{K " ] ~ ^7 ^ (A^o (448) 

which can be rewritten as a sum of eight terms: 



cn,n-i Vu^ 



h N - 1 9 3 (x N )9 3 (x N - 1 ) N(X - fl) 

y r f^ f^ Nx + m ~ ku^)6 3 {Nx + T]V + k U(X) ) (i v r , . , , , , ;i , h , , 
^-^ ' ' 9i(eu — k,u 00 )9 1 (t]v + ku^) 



e,r),K=±l 

(4.49) 
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We will see below that not all the terms contribute to the same order. 



4.4.1 Regime |A - p\ ~ 0(1/N) 

The eight terms of (4.49) can be rewritten as four combinations of the type: 
sm(iV7r(C(A) ±C(/i))) N ( X -p) (iV7r(C(A) ± ((p))) _ ^ 

In the limit |A — p\ small, i.e. \u — v\ small, the terms with a cosine will be proportional 
to derivatives of g(u, v), and there will be an overall factor. Similarly, the term with 
a sine and a + sign will be proportional to a derivative of f(u,v) and will be of order 
1/N. Only the term proportional to siniV7r J? 1 p(z)dz can balance the 1/N factor, and 
is dominant in the short range regime. After calculation we get: 

siniV7r (T p(z)dz 

K{\p) ~ Ar I K / 4.50 

v ^ J \\-nho(i/N) Nn(X-p) K ' 

As expected we have 

K(A,A)=p(A) (4.51) 
and we recover the universal short range correlation function: 

4.4.2 Long range regime, smoothed oscillations 

When | A — p\ ~ 0(1), K(X,p) has high frequency oscillations, and only a smoothed 
correlation function obtained by averaging the oscillations can be observed. 

Recall that the connected 2-point correlation function is related to K 2 by (4.9): 

p 2c (X,p) = -K(X,p) 2 (4.53) 

with K(X,p) given by 4.49. 

Smoothing out the oscillations amounts to kill all terms containing some e tN7T( > in 
the square of eq 4.49, we thus have: 



-2u'v'c 2 



~t7T\ Y2 _ C NN-1 

"^-A 2 (^)^(^-i)iv 2 (A-^) 2 £ ^ ±1 KlK2 

9 3 (Nx + eu — Kiu^O^Nx — eu — k 2 u O0 )6 3 {Nx + rjv + Kiu^O^Nx — rjv + ^u^) 
8 1 (eu - K^u^O^-eu - K 2 u oo )0 1 (r]v + Kiu 0O )6 1 (-r]v + k 2 -Uoo) 

(4.54) 
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Using that (see appendix A, and [23, 25, 24]): 

A _ 2A e l (u-v)e l (u + v)ef(2u 00 ) 

we get (eij = -eji = ±1): 

2«W 1 2 (0) 



NHn 2 9l( y Nx)9 2 1 ( y 2u c 



( O^U - Uoo)^l(« + ^oo)0l(^ - «oo)^l(u + "ocA / \ 

777 tttt — ; — s > J {^ij^ki + Ul^kj) (4.56) 

V e 1 (u-v)e 1 (u + v) J ,.^ ±i 

9s(Nx + u — iu 00 )9^{Nx — u — ku 00 )9^{Nx + v — ju 00 )9^{Nx — v — Z-Uoo) 
9i(u - iu^O^-u - ku^O^v -ju OQ )9i(-v - lu^) 

We see that 4.56 has no pole when u = ±-Uoo, it can have (double) poles only when 
u = ±v. Thus, 4.56 can be rewritten in terms of Weirstrass functions of u + v and 
u — v: 



K{\ (i) 2 = u'v' (C lP (« -v) + C 2 p{u + v)- 2S) 

Taking u = v and u = —v in 4.56, we find that the residues are C\ — C% — 1, and 
taking a particular value of u and v, we find the constant S, equal to what we had in 
3.31: 

K(\,fi) 2 = "—^2 u ' v ' + ^ u + v ^~ 2 ^ Nx + ^) ( 4 - 57 ) 

and we recover the result 3.31 found in section 3.6. 



5 Conclusions 

In this article, we have solved the puzzle raised by [11, 13] and understood why the 
naive mean- field method [10] and the orthogonal polynomial ansatz [11, 12] approach 
used in the symmetric case disagree. 

We have proven here that this effect has nothing to do with a Z2 symmetry breaking, 
as it was sometimes assumed [11], it is general as soon as the support of the density is 
not connected. 

The apparent paradox comes from the fact that when the support of eigenvalues 
is not connected, the free energy admits no large N expansion in powers of 1/N 2 
(topological expansion [9]). This means that the free energy in the multi-cut case 
is not given by a topological expansion, i.e. the sum of diagrams with a weight N x 
(x=Euler Characteristic of the diagram). 
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The explanation lies in the discreteness of the number of eigenvalues. For instance 
in the symmetric 2-cut case, the classical approach assumes that the minimum of the 
free energy is reached when one half (x = 1/2) of the eigenvalues are in each cut. 
Obviously, this minimum is never reached when the total number of eigenvalues is odd, 
and in general, the result depends on the fractional part of Nx. 

At leading order in N only, the free energy is correctly given by the classical saddle 
point limit [16, 10], but the first order in N is not sufficient to determine the 2-point 
(or higher) correlation function. 

Here we have computed explicitly the two-point connected correlation function. It 
contains a universal part depending only on the number of cuts, which was obtained 
by [10], and contains in addition, a non universal term quasiperiodic in N [22, 14]. 

Let us stress that our calculation holds for any potential, not necessarily symmetric, 
and it can also be generalized to a potential with complex coefficients (appendix. B), 
and to an arbitrary number of cuts (appendix. C). 

We have also reobtained directly the asymptotic expressions for the orthogonal 
polynomials [14], which allows in principle through the Darboux-Christoffel theorem 
(eq. 4.10) to compute any correlation function of any number of eigenvalues in the 
short or long range domain (and one can smooth it afterwards). 

The orthogonal polynomial approach may in turn be used for other random matrix 
ensembles, and it would be interesting to apply our results to orthogonal or symplectic 
ensembles [2]. 

The authors are thankful to K. Mallick for useful discussions, and to the Eurogrid 
European Network HPRN-CT-1999-00161 for supporting part of the work. They also 
thank E. Kanzieper, O. Lechtenfeld and G. Akemann for their interest and for pointing 
some missing references. 



Appendix A A few useful identities on elliptic func- 



Here we collect a few useful identities on elliptic functions used through the paper. For 
details see [23, 24, 25]. We start from 



tions 



(7(A) = (A-a)(A-6)(A-c)(A-d) 



a < b < c < d 



(A.l) 



and the map from the complex plane to the torus 




(A.2) 
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u(c)= "t/2 



u(b)=l/2+ t/2 



a b 



c d 



u(d)=0 




u(a)=l/2 



Figure 4: the upper half-plane is mapped onto a rectangle (1/2, t/2) 



where the half period K is 

" c dz 



Jb 



KT[m] = 



7r^(0|r) 



^/(c — a)(d — 6) 



(A.3) 



K[m] is the standard complete elliptic integral [23, 24, 25], with the modulus m equal 
to the biratio of the four points a, b, c, d: 



m = 



(d — a)(c — 6) 



(d-6)(c-a) 

m is related to the modular parameter r of the torus by: 



(A.4) 



m = e 



J7TT O V 2 I / 



and conversely 



.K[l-m] 

K\m] 



where we have used the Jacobi theta functions: 



eAzW) 



e 1 (z) = -i Yl 

reZ+l/2 



<f 2 e **rz with g = e « 

1 T, 



and 3 (*|r) = 3 (*) = 4 3 e OT2 0^ + - + -|r) 



(A.5) 



(A.6) 



(A.7) 



With this mapping u(X) between the A complex plane and the periodic rectangle 
of sides (1,t), we have: 



u(d) = 



u(a) = 



u{b) = 



1 + r 



u{c) = 



2 1 w 2 ' w 2 
The inverse mapping can be written in terms of theta functions: 

0;(o) 0?(u)0i(2 Woo ) 



A-d 



2X 6>i(« + Uoo)9i(u - Uoo^KUc 



g(0) 9 1 (2u)9 1 (2u OQ ) 
AK* 6\{u - UoD )6l(u + Uoo ) 

28 



u(oo) = Uoo 

(A.8) 



(A.9) 
(A.10) 



and in terms of the usual trigonometric elliptic functions sn, cn, dn [23, 24, 25] that 
we normalize to have periods 1 and r, i.e. 



sn(u) = sn(2 K[m]u\m) , dn(tt) = dn(2 K[m\u\m) , ... (A. 11' 
one has 

sn 2 (w) A — d cn 2 (u) A — a dn 2 (w) A — b 



sn 2 (-Uoo) A - c cn 2 (-Uoo) A - c dn 2 (-Uoo) A - c 



(A.12) 



; , . dn 2 (w) 



\-c=(d-c) , X-b=(d-b) (A.13) 

v ' y sn 2 (u) v ' y sn 2 (u) v ' 

sn 2 (uoo) sn 2 (« 00 ) 

X-a = (d-a) Cn2( l , X-d={d-c) d -^ (A.14) 

sn 2 (« 00 ) sn 2 (« 00 ) 

/ — ttt , , w, N ld — b m(u) cn(u) dn(-u) ,. 
y^A) = d - c)(d - a) J- — - v V (A.15) 

VC— a sn 2 (u) V 

y sri 2 ^) J 

and Moo is related to a, b, c, d by any of the following relations: 

2 / \ ^ ^ 2 / \ ^ ^ 2 / \ ^ ^ /a \ 

sn (Woo) = , cn (Woo) = , dn (u^) = - — - (A.16) 

a — a d — a d — b 

Appendix B Complex potentials 

The case of complex potentials, that is to say of a polynomial potential V(A) with 
complex coefficients, is interesting for some applications of the matrix models to 2 
dimensional gravity and when studying their connections with integrable hierarchies. 
In this case, the mean field large N solution is known to be given by a continuous 
distribution of the eigenvalues along arcs in the complex plane [26]. 

In this appendix we show that our results are only slightly modified in this case. 

In the two-cut case, we can repeat the analysis of sect. 2. We fix x = n\/N (the 
proportion of e.v. in the first cut). The resolvent is still of the form 2.28, with the 
polynomial M and the end points a, b, c, d determined by the constraints 2. 18 and 2.29, 
but they are no more real in general, as well as the resulting mean-field free energy 
F (x). 

If we now repeat the calculation of sect. 2. 2. 3 we cannot use a saddle-point approx- 
imation for the sum over n\ by expanding F (x) around the saddle point x which is 
the true extremum of F . 

a>> = < B1 > 
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Indeed this extremum is at a finite non-zero distance of the real axis, i.e. Im(xo) = 
0(1), while the method of sect. 2. 2. 3 is valid only if Im(a;o) = 0(1/ N). However, since 
iV is integer, we can expand F around any x k , provided that 

F^x k ) = 2m ^ , keZ (B.2) 

since the dangerous oscillating term e~ N2 ^ x ~ Xk ^ F o^ i s then a constant for x = n/N, 
neZ. 

Therefore, as in [26] we have to consider the real pseudo saddle-point x c such that 
Re(i^(a; c )) = with lm(x c ) = (B.3) 

and denote 

Ac = i F ° (Xc) = ^ Im ( F o(^c)) (B.4) 
We expand F around some x k defined by Eq. (B.2) and such that 

x k -x c = 0(N~ l ) (B.5) 

and we get for the total free energy (by exactly the same calculation as in sect. 2. 2. 3) 
F = N 2 F (x k ) - \n(e 3 (Nx k )) 

(B.6) 

+ F 1 (x c ) + \ ln(27r^'(x c )) + 0(N~ 2 ) 



where 6*3 is the theta function with modular parameter 

2i7T 



(B.7) 



F{(x c ) 

Only the first two terms are important for calculating the two-point functions and the 
orthogonal polynomials in the large iV limit. This leading term does not depend on k. 
Indeed we have 

F^(x k ) - F^(x c ) = (x k - x c )F£(x c ) + 0(N~ l ) (B.8) 

hence 

(x k - x c ) = 1 r(k - NA C ) + 0(N- 2 ) (B.9) 

and using the periodicity relations of 9 3 we can rewrite the leading term for the free 
energy as 

17T 

N 2 F (x k ) - \n(9 3 (Nx k )) = N 2 F (n c /N) u\ - ln(^K)) (B.10) 

T 

with 

n c = E[Nx c ] , u c = [Nx c }-t[NA c ] (B.ll) 
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where E[w] is the integer part of u (largest integer smaller than u) and [u] — u — E[u] 
is the fractional part of u. This does not depend on k up to negligible terms of order 
0{N~ 1 ) (provided that condition B.5 for k holds). 

One can now repeat the calculation of sect. 3 for the 2-point function. Nothing is 
changed but we simply have to replace x c by X}~ in the intermediate steps and to use 
Eq. (B.9) at the end of the calculation. This amounts to replace the x c in the elliptic 
function sn 2 by x c — rA c . The final result for the two-point resolvent is 



w c (A,/i) 



1 



4(A-/i) 



2 



' " " ^-a)^-b)(X-c)(X-d) 1 " (A ~ /0 



{c-a){d-b) 2 



4 V ^(A) V ^) 



sn (jV(x c — rA c ) + |) (B.12) 



Similar results holds for the orthogonal polynomials. We simply have to consider the 
end-points a, b, c, d for the mean-field real parameter x c and to make the replacement 

Nx c -> N(x c + tA c ) (B.13) 

in the elliptic functions involving Nx c . In any case these terms depend only on the 
fractional parts of Nx c and of iVA c . 

A final interesting remark on the periodicity properties of the non-universal term 
sn 2 (iV(£ c + tA c )) can be made. From the definition B.4 A c corresponds to a "phase 
shift" between the two arcs where the density of e.v. is non zero. 

where NT a is the (constant) effective potential on the arc a. The two periods 1 and 
r of the sn 2 function correspond respectively in term of eigenvalues to (i) transfer a 
single e.v. from the first arc to the second one (5Nx = ±1), (ii) or to shift the phase 
between the two arcs by 2-7T (SNA C = ±1). 

Appendix C Multicut Case 

Consider now a support of eigenvalues split into s intervals: 

C = CiU ...DC s (C.l) 

Let rii be the number of eigenvalues in each Ci, and Xi = rii/N the occupation ratio, 
which we denote collectively as a vector: 

Xi= dXp(X) , x = (x 1 ,...,x a -i) (C.2) 
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Note that only s — 1 of them are independent since X\ + . . . + x s — 1. 

As in the two-cut case (eq. 2.36), the free energy at fixed n admits a topological 
large N expansion: 



F[V; n] = N 2 F [V; x] + N°Fi [V; x] + 0{1/N 2 ) 
and as in 2.38, the partition function can be written as a sum over n: 

Z = e~ F = J2 e~ F[VA 

n 

The sum is dominated by the vicinity of the extremum x c of F [V; x\: 
Z~J2 e-^ ^^- 3 ^' 1 ^-^) where 

n 

risthes — lxs — 1 matrix defined by: 

-i 1 d 2 F 



(C.3) 



(C.4) 



= (C.5) 



13 2m dxidxj 



Then, the summation over n yields: 

Z ~ e -^o[v,s c ] e{NSc \ T) 

where 6{u\r) is Riemann's theta function [27] in genus s — 1: 

6{u\t) = 6{u) ^^ e -i<n-u).r-\n-u) = ^^n.rn^ 



2inn.il 



(C.6) 



(C.7) 



(C.8) 



where r is as — lxs — 1 matrix, u is a s — 1 component vector, and n is a vector with 
integer coordinates. 

The 9 function obeys the relations (k being an arbitrary integer vector): 



9{u + k) 



u 



6(u + rk) = e -™W-k+k.rk) 9 ( 



u 



d(-u) = 6{u) (C.9) 



Eventually the free energy at leading orders in iV is: 



F ~ N 2 



F [V, x c \ - ± ^0(Nx c \t) + j^F^xJ + ... 



(CIO) 



It is now straightforward but lengthy to rederive the 2-point correlation function 
and the orthogonal polynomials from C.10. One needs to differentiate C.10 with respect 
to variations of the potential as in 3.3 or 4.16, and express the hyperelliptical functions 
involved in the calculation through prime forms (hyperelliptical generalization of the 
6*i function) [27]. One should thus obtain expressions similar to those of [14]. 
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